Spectral correlations of individual quantum graphs 
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I. INTRODUCTION 

The spectral fluctuations of individual complex 
(chaotic) quantum systems are universal and can be 
described in terms of Wigner-Dyson random matrix 
theoryd Q. (See also H S H B 

and references 

therein). For classically chaotic systems this empirical 
statement was promoted to a conjecture by Bohigas, Gi- 
annoni and Schmit 0. (See also While, however, 

there is enormous experimental and numerical evidence 
in support of this conjecture the physical basis of 
universality is not yet fully understood theoretically. 

To date, the most advanced approach in develop- 
ing correspondences between spectral statistic and non- 
linear dynamics is semiclassical analysis. Beginning with 
Berry's seminal work 12j it became understood that in- 
formation on spectral correlations is stored in action cor- 
relations of classical periodic orbits (see also 0|.) Going 
beyond the 'diagonal' approximation [l2j wherein only 
identical (and mutually time-reversed) orbits contribut- 
ing to the Gutzwiller double sum are taken into ac- 
count, a hierarchy of ever more complex expansions in 
orbit pairs has been constructed [13,13 US- In this way, 
it was shown that to all orders in an expansion in the 
ratio T = t/tn the short time {t < tn) behavior of the 
spectral form factor Kir) of imiformly hyperbolic quan- 
tum systems coincides with the universal predictions of 
random matrix theory (RMT). (Here, tn — ^§ denotes 
the Heisenberg time and AE is the mean level spacing.) 
However, in view of the fact that at r = 1, the function 
K(t) contains an essential singularity, it is presently not 
clear how to extend this expansion to times larger than 
the Heisenberg time. 

Some time ago, a field theoretical approach to quan- 
tum chaos — dubbed the ballistic cr-model — has been 
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introduced |l8l | as an alternative to semiclassical expan- 
sions. The most promising aspect of this development 
is that in field theory the full information on universal 
RMT correlations is obtained in a very simple manner, 
viz. by integration over globally uniform 'mean field' 
configurations; universality of chaotic quantum systems 
is proven, once it has been shown that at sufficiently low 
energies (long times) fluctuations become negligible and 
the field theory indeed reduces to its mean field sector. 
Unfortunately, however, it has so far not been possible to 
demonstrate this reduction in a truly convincing manner. 
(The situation is much better in the field of disordered 
chaotic systems: It has been known for some time that 
at low energies disordered systems exhibit RMT spectral 
correlations upon configurational averaging. This type 
of universality has been proven 7] by field theoretical 
methods similar to those mentioned above.) 

Motivated by the lack of universality proofs for generic 
quantum systems with underlying Hamiltonian chaos, we 
have recently considered the spectral properties of quan- 
tum graphs |l9ll. (For the general theory of quantum 
graphs, see |23| and references therein.) Quantum graphs 
differ from generic Hamiltonian systems in two crucial as- 
pects: First, the classical dynamics on the graph is not 
deterministic. It is rather described by a Markov pro- 
cess. Second, quantum graphs are 'semiclassically exact' 
in that their spectrum can be exactly described in terms 
of trace formulae. In spite of these differences, quantum 
graphs display much of the behavior of generic hyperbolic 
quantum systems (while being not quite as defiant 
to analytical treatment than these.) 

Earlier work on universal spectral statistics in quantum 
graphs was based on periodic orbit summation schemes 
similar in spirit to the semiclassical approach to Hamil- 
tonian systems. Specifically, Berkolaiko et al. |22J de- 
veloped a perturbative diagrammatic language to ana- 
lyze the periodic-orbit expansions of spectral correla- 
tion functions beyond the diagonal approximation. Tan- 
ner 01 analyzed the structure of the semiclassical ex- 
pansion to conjecture criteria for the presence of univer- 



2 



sal correlations on graphs. He also established connec- 
tions between universality and the decay rates of clas- 
sical Markovian dynamics of the system (for details see 
appendix . 

While all building blocks of semiclassical analysis on 
graphs are known l24j . and a complete summation 
over all orbit pairs may be in reach, semiclassics on 
graphs is subject to the same limitations as in Hamilto- 
nian systems. In particular, it is not clear how to extend 
its domain of applicability to times beyond the Heisen- 
berg time. In view of these difficulties, we have developed 
an alternative approach 19] which is based on field theo- 
retical methods and avoids diagrammatic resummations 
altogether. Rather, it is based on two alternative pieces 
of input, both of which have been discussed separately 
before: 

i. The exact equivalence of a spectral average for a 
quantum graph with incommensurate bond lengths 
to an av erag e over a certain ensemble of unitary 
matrices |2ll I23ll25| . 

ii. An exact mapping of the phase-averaged spec- 
tral correlation functions onto a variant of the su- 
persymmetric cr-model by an integral transform 
known as the color-flavor transformation . 

The synthesis of i. and ii. leads to a formulation simi- 
lar in spirit to the 'ballistic cr-model' yet not burdened by 
the technical problems of that approach. It is the purpose 
of this paper to give a detailed account of this theory, and 
to discuss a number of generalizations. Specifically, we 
will discuss the crossover between systems of conserved 
(orthogonal symmetry) and broken (unitary symmetry) 
time-reversal invariance, and we will consider the case of 
broken spin rotation invariance (symplectic symmetry.) 

The paper is organized as follows: In Section ^1 we 
give a short introduction to quantum graphs. We discuss 
the relevant quantization conditions, spectral correlators, 
and the meaning of incommensurate bond lengths. The 
supersymmetry approach to the spectral two-point cor- 
relation function is discussed in Section ITTll In Section 
II VI we subject the supersymmetric generating functional 
to a stationary phase analysis. We show under which 
conditions the field theory can be reduced to a 'mean 
field' theory of RMT-type correlations. We also discuss 
the crossover between graphs of orthogonal and unitary 
symmetry. Quantum graphs belonging to the symplec- 
tic symmetry class are discussed in appendix |nj and an 
outline of the proof of the color-flavor transformation is 
given in appendix lUl 

II. QUANTUM GRAPHS 

A. Generalities 

A finite graph Q consists of V vertices which are con- 
nected by B bonds. The V x V connectivity matrix is 



defined by 

Cij — # {bonds connecting the vertices i and j} . (1) 

A graph is simple if for all i, j, Cij £ [0, 1] (no parallel 
connections) and Ci^i — (no loops). The number of 
bonds is B = ^ "^Jj^i The valency of a vertex i 

is the number of bonds connected to it Vi — X]]=i ^i-j- 
A graph is called 'connected' if it cannot be split into 
disjoint sub-graphs. With only slight loss of general- 
ity [23|, we will focus on the case of simple connected 
graphs throughout (cf. Fig. ^for a schematic.) 




FIG. 1: Schematic of a fraction of a (planar) simple quantum 
graph and of the notation used in the text. 

We denote a bond connecting the vertices i and j by 
b = [ij). The notation {ij) and the letter b will be used 
whenever we refer to bonds without specifying a direc- 
tion: b e {ij) = [ji). A directed bond /3 = (6, c?) consists 
of a bond b and a direction on the bond which will be 
denoted by a direction index d = ±1. For b = (ij) and 
i < j we set d — +1 for the direction i —f j and d = —1 
on the opposite direction. 

The position a; of a point on the graph is determined 
by specifying its bond 6, and its distance Xf, G [0, L\,] from 
the adjacent vertex with the smaller index. The length 
of a bond is denoted by Lf,. Throughout, we will assume 
the bond lengths to be incommensurable (or rationally 
independent) in the sense that there is no non-vanishing 
set of integers mf, G K such that "^fc^fc = 0. 

The Schrodinger operator on Q is defined by one- 
dimensional Laplacians on the bonds, and a set of ver- 
tex boundary conditions establishing self-adjointness. Its 
wave functions ^!{x) are complex valued, piecewise con- 
tinuous and bounded functions. Writing "^{x) = ipb^Xf,) 
for X = Xb d (0, Lf,), the solutions of the stationary 
Schrodinger equation at a given wave number fc > have 
the form 
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where ab^d=±i are the complex ampHtudes of 'right' 
{d — +1) and 'left' {d — —I) propagating waves on the 
bond, and Af, are constant phases generated by optional 
magnetic fluxes threading the plaquettes of the graph. 
To characterize the vertex boundary conditions, we in- 
troduce (fc-independent) Vi x Vi vertex scattering matri- 
ces connecting incoming waves to outgoing waves 
at i. (6 and b' run over bonds connected to i). Denoting 
the outgoing/incoming direction on bond b by dout/c^in, 
these matrices are defined by the equation 



"Mo 



b 



(3) 



To represent this equation in a more concise form, we 
combine all amplitudes abd into a 2i?-dimensional vector 
a. In this notation, 

d^UB{k)a (4) 
where the 2B x 2B quantum map is given by 

Ueik) ^T{k)ST{k), (5) 



the diagonal matrix T{k) 



bd.b' d' 



bd,b' d'G 



i{k+dAb)Lb/2 



describes the propagation along half a bond, and 



bd,b'd' 







ib',d') 

else, 



(6) 



combines all boundary conditions at the vertices into a 
single scattering matrix. The equivalent of the quantum 
map in a Hamiltonian system is a quantized Poincare 
map. 

The boundary condition can be fulfilled only for 
discrete set of wave numbers fc„. These numbers define 
the spectrum of the quantum graph. For fc > the quan- 
tization condition Q is equivalent to the vanishing of the 
spectral determinant 



£,{k) EE det{l^UB{k)) . 



(7) 



Thus ^(fc) = if and only if fc = fc„ > is in the spec- 
trum. 

By way of example, we mention two frequently em- 
ployed families of boundary conditions: so-called Neu- 
mann boundary conditions 29] correspond to 



Sb.b'Vi 



(8) 



For large valencies Vi the non-diagonal terms are much 
smaller than the diagonal, and the back scattering term 
dominates. While on general graphs wave functions need 
not be continuous across the vertices, they are so on 
Neumann graphs . Another interesting set of bound- 
ary conditions is implemented through discrete Fourier 
transform (DFT) matrices 



^(i),DFT _ 1 



"6,6' 



exp 



2TTi 



.7r(6)7r(6') 



(9) 



where b i— > 7r(6) maps the Vi bonds connected to vertex i 
one-to-one onto the numbers 1, 2 . . . , w^. These boundary 
conditions do not imply continuity at the vertices; incom- 
ing wave packets are scattered into the outgoing bonds 
with equal probability. 



B. Time— reversal invariance 

As with Hamiltonian chaotic symmetries, quantum 
graphs of different symmetries may be identified. Specif- 
ically, quantum graphs carrying spin degrees of freedom 
(and spin-rotation invariance breaking vertex scattering 
matrices) fall into the symplectic or unitary symmetry 
class depending on whether time reversal invariance is 
broken or not. These cases will be discussed in ap- 
pendix^ In the absence of spin, we need to distinguish 
between graphs with broken (unitary symmetry or sym- 
metry class A in the notation of jS^) or conserved (or- 
thogonal symmetry or symmetry class AT) time reversal 
invariance. 

A quantum system is time-reversal invariant if its 
Hamiltonian H commutes with an anti-unitary time- 
reversal operator T, [H,T] = 0. For spinless systems, 
T is an involutory operator, — 15]. The condition 
[H, T] — restricts the form of both the bond propa- 
gation matrix T(k) and the vertex scattering matrices 
S(fc). In non-time reversal invariant systems, these ma- 
trices obey no conditions other than unitarity. However, 
for conserved time reversal invariance, and with a defini- 
tive choice of the time-reversal operator T, all vertex 
scattering matrices have to be symmetric E*^'^ = S'^'^ , 
i.e. 



(10) 



where = (Jo) Pauli matrix in direction in- 

dices d. Additionally, all magnetic phases A^ must van- 
ish. The crossover between orthogonal and unitary sym- 
metry will be discussed in IIV CI where we explore the 
consequences of a gradual switching on of magnetic phase 
factors. 



C. The density of states and spectral correlation 
functions of quantum graphs 

The density of states (DoS) of a quantum graph is 
defined as 

rf(fc) = I] '^(^ - ^«) = ^ + M(fc) (11) 



where the sum runs over the spectrum fc„ . We have writ- 
ten the DoS as a sum over a smooth part 1/ A where A 
is the mean level spacing and fluctuations 5d{k). Both 
parts allow for an explicit representation in terms of the 
quantum evolution map. For the mean (or Weyl) part 
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one obtains 
1 

A 



1 d 

27ri dk 



Indeti-Usik)) 



BL 

TT ' 



(12) 



where L = ^ Li, is the mean bond length. Note that 
the mean level spacing is constant. The fluctuations can 
be expressed through the spectral determinant pH 



5d{k) = -ilm4-ln^(fc+), 
TT dk 



(13) 



where k^ = k + ie and the limit e — > is implied. Using 
that In det A = txhi A and expanding the logarithm one 
obtains an exact Gutzwiller type trace formula 



5d{k) = -Im4 V -trZi|^(fc+) 
^ ' TT dk n 

n—l 



(14) 



expressing the DoS in terms of a sum over periodic orbits 
(periodic sequences of directed bonds.) 

We aim to explore the statistical properties of the fluc- 
tuating part of the DoS. The iV-point DoS correlation 
function is defined by an average over the complete spec- 
trum 



i?7v(si, • ■ • 
{Sd{k - 

where 



SN-l) = 
- SAT-iA) , 



.d{k + siA)Sd{k)),^, 



lim — 

K^oo K 



K 



dkf{k). 



(15) 



(16) 



Throughout, we will focus attention on the two-point 
correlation function i?2(s). The two-point function can 
be conveniently expressed as a derivative of quotients of 
spectral determinants: 



i?2(s) 

where 



1 



d2 



,,=oI^c (C(j+,J-;s))fc, (17) 



C(i+, j-;s) 



C(fc++p+f) /e(fc++P-f) 



^(fc^ 



p-h) 



(18) 



and p±b - (±s/2 - j±) A, p±f = (±.s/2 + i±) A. (Higher 
order correlation functions may be obtained in a similar 
manner from generating functions involving additional 
quotients of spectral determinants.) 

For later reference, we recall that the RMT two-point 
correlation functions are given by 



i?«u^(.) =5{s) 



■ 2 

sm ITS 



R^OB^S) =i?2^™(.) 



2 V^y"^ (19) 

(7r|s| cosTTS — sin7r|s|) (2Si(7r|s|) — tt) 
2^ ' 



where Si (a;) = /^^ dx' '^'"f' is the sine integral. 



We also notice that the statistical properties of the 
graph may be characterized by correlation functions dif- 
ferent yet closely allied to the correlation functions in- 
troduced above: for any value of k the quantum map Ub 
possess a set of 2B 'eigenphases' e^**'^*^^ {I — 1, ... , 2B) 
on the unit-circle. At fixed k the density of phases is 
given by 



2B 



Pk{e) ^Y.^2-n{e - ei{k)) 



where 52tt{&) is the 27r-periodic delta-function. The sta- 
tistical properties of this quantity are defined by averag- 
ing over both and k. Occasionally — e.g. within the 
context of the periodic orbit approach to graphs — it is 
sometimes advantageous to consider the correlation func- 
tions of the eigenphases instead of the spectral correlators 
introduced above. Under mild conditions (weak fluctua- 
tions in the bond lengths) both types of correlators are 
equivalent in the limit _B — > oo of large graphs. In the 
following we will keep our discussion focused on the spec- 
tral correlators. With small and straight forward changes 
our theory can be applied to the eigenphase-correlations 
as well. 



D. Consequences of incommensurability 

The quantum map Usik) = T{k)ST{k) depends 
on the wave number k via the B diagonal elements 
gifcLi,/2_ Defining (l>b{k) = fcLb/2 we have a map k i-^ 
(g20i(fc)^ _ g«0B(fc)) of the wavenumber into a B-torus 
C <C^ . This map may be interpreted as a 'Hamil- 
tonian flow' where k plays the role of 'time'. As we as- 
sume incommensurable bond lengths if,, the image of the 
phase map covers the torus densely, i.e. the Hamiltonian 
flow is 'ergodic'. This in turn means that long time aver- 
ages (/c-averages) may be traded for phase space averages 
(averages over the torus or, equrvalently, independent av- 
erages over the B phases 4>b UM)' 



■{2n) 



(20) 



It is this equivalence which makes the analytical calcu- 
lation of spectral correlation functions a feasible task. 
Upon replacing ( )k { )ct>, the one-parameter fam- 
ily of matrices T{k) — > T{{(j)i,}) becomes an 'ensemble' 
of random matrices. There is a well developed analytical 
machinery designed to perform random phase averages of 
this kind. Below, we will apply the formalism of super- 
symmetry to compute the random phase averaged spec- 
tral correlation functions of the graph which, by virtue of 
the equivalence above, are strictly equivalent to the wave 
number averaged correlation functions. 
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III. NONLINEAR a MODEL FOR QUANTUM 
GRAPHS 

Consider the representation (|17(l of the two-point cor- 
relation function in terms of a double derivative of the 
quotient ( of spectral determinants. Replacing the k- 
average by a random phase average, {Qk (C)^: it is 
the purpose of the present section to derive a cr-model 
representation of the two-point correlation function. 



A. The generating function as a Gaussian 
superintegral 



Defining the supervectors 
/.si\ 



SN 

Xi 

\XN/ 



and V — (^1 



XI ...Xn), (21) 



where Si are complex commuting variables while Xi ^-nd 
Xi are independent anti-commuting numbers, the quo- 
tient of determinants of an TV x matrix Af and a (pos- 
itive) Nx N matrix A^, can be represented as a Gaussian 
integral 

= sdet A = / d{iP,^)e-^^'^. (22) 



det Ab 



Here, 



A^ 



Ab 
At 



(23) 



is a block-matrix in boson-fermion space (the two com- 
ponent space introduced by the s/x grading of ijj) and 
the measure is given by 



N 



diij), -0) = TT ^ dRe{si)dlui{s^)dxidxi 



(24) 



where / dxi Xi = 1 ^ind / dxi = 0. We wish to apply 
this relation to represent the spectral determinants (0 
in terms of Gaussian integrals. In view of our applica- 
tions below, it will be convenient to double the matrix 
dimensions j^J using 



C(fc+p) =det {UBip))det 
which leads to 

CU+J-;s) 

where 



1 T{k) 
T{k) UsipY 

fi(V;,V)(sdet(T+Tl))2e-^['^''''l 



V'- 



1 r(/c)t 

T(fc)t T^ST^ 



(25) 



(26) 



(27) 



V'- 



Here, -0 = {il'a,s,x,d,b} is a 16i?-dimensional supervector 
where, a = ± distinguishes between the retarded and the 
advanced sector of the theory (components coupling to ^ 
or , respectively) . The index s = f , b refers to complex 
commuting and anti-commuting components (determi- 
nants in the denominator and numerator, respectively), 
and a; = 1, 2 to the internal structure of the matrix kernel 
appearing in ()27|l . The matrices 



T4 



- [ 





T{p±,f] 



(28) 



are diagonal matrices in superspace containing the ap- 
propriate bond matrices T in the boson-boson/fermion- 
fermion sector. 

To account for the (optional) time-reversal invariance 
of the scattering matrix, we introduce the doublets 



V2 V'^f'V'^ 



'V2 



(■0, 



(29) 



where cr^^ 



) is the Pauli matrix in superspace. 
Notice that the lower components of 5* emanate from the 
upper component by a time reversal operations (trans- 
position followed by inversion in directional space.) For 
later reference, we note that new fields depend on each 
other through the generalized transposition 



* = r. 



The explicit definition of the matrix r is given by 



r =(Ti To 



TO =Ebal'' - ^^fa*^ 



(30) 



(31) 



where cr*'" are Pauli matrices in the newly introduced 
'time- reversal' space and -Eb/f are the projectors on the 
bosonic/fermionic sectors. However, all we will need to 
know to proceed is that r obeys the conditions 



-1 

T — T 



and 



..BF 



(32) 



The appearance of the matrix r in conjunction with a 
transposition operation suggests to introduce the gener- 
alized matrix transposition 



A^ = tA' t 



T_-l 



(33) 



Using Eq. (Eg and that [A^f = af^Aaf^, one finds 
that the generalized transposition in involutory, 

{A^y = A. (34) 

For later reference we also note that 

*+.4*_ = ^'^af^A'^^l = vi>_A^*+. (35) 
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With all these definitions, the action (|27|l now takes the 
form 



lT(fc)t S- ' 



(36) 



where the matrix structure is again in the auxiliary index 
X and we have introduced the matrices 

s± - r±5r± EE T±s%±) 

Here the matrix structure is in time-reversal space and 
the time-reversed scattering matrix has been defined 
in fTn|. 



B. The color-flavor transformation 

We are now in a position to subject the generating 
functional to the spectral average, (C)fc — {C)<t>- As 
discussed in section ITl Dl we replace T(fc) T{{(j)h}), 
whereupon the average is given by 

(C(j+,j-;s))0= J d(vi,V)sdct(r+r!)2x 

B r^A (38) 



-,-So 



b=l 



Here, 



is the phase-independent part of the action and 

Si.b - 2^+s,be'^'-^+^2.b + 2*_.2,6e-'^^*_4,f,. (40) 

So far, we have not achieved much other than represent- 
ing the spectral determinants by a complicated Gaussian 
integral, averaged over phase degrees of freedom. The 
most important step in our analysis will now be to subject 
the generating function to an integral transform known 
as the color-flavor transformation |2g . The color-flavor 
transformation amounts to a replacement of the phase- 
integral by an integral over a new degree of freedom, Z . 
Much better than the original degrees of freedom, the 
Z-field will be suited to describe the low energy physics 
of the system. 

In a variant adopted to the present context (a single 
'color' and F 'flavors') the color-flavor transformation 
assumes the form 



2-n 



j d{Z, Z)sdet (l - ZZ^ ^vlzn^+u^zu^ ^ ^^^^ 



where r]± and i/± are arbitrary 2F dimensional supervec- 
tors and Z, Z are 2F-dimensional supermatrices. The 
boson-boson and fermion-fermion block of these super- 
matrices are related by ^bb = -^bb: -^ff = "■^ff i while the 
entries of the fermion-boson and boson-fermion blocks 
are independent anti-commuting integration variables. 
The integration d{Z , Z) runs over all independent matrix 
elements of Z and Z such that all eigenvalues of Z\^\:,zl^^ 
are less than unity and the measure is normalized such 
that 



/ 



(i(Z, Z) sdet(l - ZZ) = 1. 



(42) 



51 T 
— ► — »- 



z 



pt 



FIG. 2; On the physical interpretation of the color-flavor 
transformation. Explanation, see text. 

We apply the color-flavor transformation B times - 
once for each phase 0b. As a result, we obtain a i3-fold 
integral over supermatrices Zb. There are four flavors 
(direction index d — ±1 and time-reversal index t = 1,2). 
We combine all matrices Zt, (Zt,) into a single block- 
diagonal 8i?-dimensional supermatrix Z (Z) such that 



Zbdts,b'd't's' — SbM'Zbdts,d't's'- (43) 

The averaged generating function now has the form 

{aj+,J-;s)) -sdet(T+T!)2 [ d(^,^) / d{Z,Z) 

J J (44) 

sdct(l - ZZ) e-s(*.*.^.^) 

where 



S{^,^,Z,Z)=^, + 



(45) 



and we used 2^xZ^i = ^iZ^i + ^iZ^^fi, 2^22^2 = 
^2^*2 + 5'2^'^*2- Here, the indices 1, 2 refer to the 
auxiliary index x, and the matrix structure is in ad- 
vanced/retarded space. Integrating the Gaussian fields 
^' and we arrive at the (exact) representation 



(C(j+,j-;s))= / d(Z,Z)e-^(^'^) 



(46) 
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where the action is given by 



S{Z, Z)^~ strlog (1 - ZZ) + -strlog (1 - Z^ Z) 



1 



strlog (1- 5+ 



(47) 



(Note, that the prefactor sdet (T+tI)^ has canceled out.) 

Before carrying on, let us pause to discuss the ad- 
vantage gained by switching to the Z-representation. 
Consider the bilinears i/'(+/-),s,x,d,6e(+/~)"^'"i/'(+/-),s,x,d,b 
appearing as building blocks of the original phase- 
representation. Loosely identifying i'+z- as re- 
tarded/advanced wave function amplitudes, these prod- 
ucts describe the scattering of single particle states off 
phase fluctuations. Due to the effective randomness of 
the phases they fluctuate in a wild and non-controllable 
manner (see Fig. |21 for a cartoon of the propagation of 
a retarded [upper line] and advanced [lower line] wave 
function in space: a rapid succession of scattering events 
[the vertical dashed lines] leads to strong fluctuations.) 
Technically, this means that the original representation 
defies controlled evaluation schemes (such as mean field 
approximations and the like.) 

In contradistinction, the Z-field enters the theory as 
~ ^ +,s,tA,bZ},^ss' ,tf ,dd''^ -,s' ,t' ,d' ,b, i-e. through structures 
that couple retarded and advanced amplitudes locally in 
space. While (prior to the phase averaging) each of the 
^+/- amplitudes individually was a rapidly fluctuating 
contribution, the product ^_|_^_ contains benign, slowly 
fluctuating contributions. This is because the phase 
exp(i(/<) picked up by the retarded amplitude may can- 
cel against the phase exp(— i(/)) carried by the advanced 
amplitude. In a semiclassical manner of speaking, this 
happens if the two amplitudes propagate along Feynman 
paths locally correlated in space. The advantage of the 
Z-representation is that it selects precisely these slowly 
fluctuating, spatially correlated bilinears which survive 
the averaging over phases. In Fig. |21 the Z-fields are 
indicated by vertical ovals. Wave function amplitudes 
qualifying to form a slowly fluctuating couple may carry 
different time-reversal and directional indices which ex- 
plains the mafra-structure of Z in these index spaces. At 
any rate, the structure of the color-flavor transformed 
theory indicates that the Z-integral will be compara- 
tively benign and amenable to stationary phase treat- 
ment. 



IV. SADDLE POINT ANALYSIS AND 
UNIVERSALITY 

The action provides for an exact representation of 
the generating functional of an individual graph. While 
the integral over Z cannot be done in closed form, it turns 
out to be ideally suited to a mean field treatment. In the 
following, we will formulate the mean field analysis and 



explore under which conditions the theory reduces to one 
that predicts universal GOE statistics. (We assume time 
reversal invariance throughout.) 

Our strategy will be to first identify uniform 'zero 
mode' solutions to the mean field equations, and the cor- 
responding mean field action. We will find that the in- 
tegral over the reduced action generates an exact RMT 
expression for the spectral determinants. In a second 
step we proceed to investigate the validity of the zero 
mode approximation, i.e. we will explore under which 
conditions corrections to the RMT result vanish in the 
semiclassical limit B —t oo. 



A. The zero-mode and universality 

We begin by expanding the full action to linear order 
in the sources p± — A(±s — <j^^ j±) 



S{Z,Z) 



--So{Z,Z) 

P+{CS 



I 

-str- 

4 



+ TStr- 

4 



SOfZ-'S^Z 



1 - 5Z^5tZ 



(48) 



1 - S^ZSZ^ 



0{{pLf 



Here, Sq{Z, Z) is obtained from (|47|l by replacing S± 
S, and Cbb' = Sbb'Lf, contains the bond length on its diag- 
onal. Since we are only interested in spectral fluctuations 
on the scale of the mean level spacing p± ^ 0{B^^) im- 
plying that higher orders in the expansion in p± vanish 
in the limit B ^ oo. At this point we have to assume 
moderate bond length fluctuations such that Lb/L <C B. 

To identify the mean field configurations of the theory, 
we differentiate the action Sq w.r.t. Z and obtain 



1 



1 - ZZ 



-Z 



1 



1 - Z-'Z 







This equation is solved by 



Z = Z^ 



(49) 



(50) 



Differentiating w.r.t. Z and using H5()|l a second saddle- 
point equation assumes the form 



1 



1-ZZ 2l-SZS^Z 
1 1 



-^SZS''- 



(51) 



2 1 -5^Z5t^Z 



-^S'^ZS''^ = 0. 



This equation is solved by all field configurations that 
commute with the scattering operators, i.e. 



Zab,dtsA't' s' —Sdd'Yts,t's' 
Zab,dtsA't' s' —Sdd'Yts,t's' , 



(52) 



which corresponds to equidistribution on the set of di- 
rected bonds. The symmetry condition Z = Z'^ obtained 
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from the first saddle-point equation implies Y = Y"^ 
where = tqY^Tq^ and the matrix tq has been defined 
in The commuting parts of these matrices obey 

^bb = ^H-d Yff = —Y^ while the non-commuting 
entries are all independent integration variables. The 
fermion-fermion part is integrated over IR4 ~ while 
boson-boson part is restricted to the compact region 
where all eigenvalues of l^^b^bb are less than unity. 

Reducing the action H48I) to the zero-mode the first 
contribution vanishes exactly So{Zo, Zq) — while the 
remaining term becomes 



For the purposes of our discussion, it will be sufficient 
to consider the expansion of the exact action H48|) to sec- 
ond order in the fields Z , 

s'^^HzJ) =si^\z, z) 



stv {p+{CS + SC)Z-S^z) (58) 



-str(p_(£5t , 



where 



^GOE(^ y) = +^^str^±I^ - ^^str^^, (53) 
' ' A l^YY A \-YY 



s'i\z,Z) = sli [ ZZ -^Z^Z - ^SIZS+Z^ ) . (59) 



Restricting the integration to the zero mode sector, we 
obtain 



(C(j+,j-;s))=^C"^(j+,j-;s)> 



d{Y,Y) 



(54) 

where the denotation C'^ce indicates that the matrix in- 
tegral over Y obtains but an exact representation of the 
GOE correlation function. To represent the integral on 
the r.h.s. in a more widely recognizable form, let us de- 
fine the 8x8 supermatrix 



Q 



1 Y 
Y 1 



1 Y 
Y 1 



1 + 2YY/{1 - YY) -2^7(1 - YY) 
2Y/{1-YY) -\-2YYl{\~YY\ 



(55) 



where = ( q _ j j . It is then a straightforward mat- 
ter to show that the action S{Y , Y) takes the form of 
Efetov's action .7| for the GOE correlation function 



C°oi=(j+,j_;s) = y dQe 



tS{Q) 



where the measure is given by dQ = d{Y , Y), 



S{Q) = -str(g-I],)e 



(56) 



(57) 



and e = (^^q p_^- For a discussion of the integral 

(|56|l . and the ways random matrix predictions are ob- 
tained by integration over Q, we refer to the textbook Q. 



B. Validity of the saddle— point approximation 



Physically, the quadratic action describes the joint prop- 
agation of a retarded and an advanced Feynman ampli- 
tude along the same path in configuration space. (This 
is a generic feature of second order expansions to non- 
linear cr-models of disordered and chaotic systems. For 
a discussion of this point, we refer to Ref. Q.) It thus 
carries information similar to that obtained from the di- 
agonal approximation to semiclassics. The second order 
expansion is justified if the fluctuations of the flelds Z 
are massively damped (in the sense that the matrix el- 
ements of Z effectively contributing to the integral are 
much smaller than unity.) Under these conditions, the 
integration over matrix elements of Z may be extended 
to infinity and we obtain a genuine Gaussian integral. 

The eigenvalues of the quadratic form appearing in 
S'*'^'' at p± =0 determine the damping mj — or the 
mass, in a field theoretical jargon — inhibiting fluctua- 
tions of the eigenmodes Zj. As indicated by its name, the 
zero-mode Zq carries zero mass. Within the quadratic 
approximation, the correlation function assumes the form 
(C(J+,j-;s)) ~ (C(J+,j-;s))(2) =n,/„with the Gaus- 
sian integrals 



JdiZ„Z,)eM-sir{m,Z,Z, 



+ i{ni — l)L{p+ZjZj — p^ZjZj)}] 

(60) 



where Zj, and Zj are 4x4 supermatrices obeying the 

ubiquitous condition Zjbb = ^j*bb ^^'^ ^jff = 
also assumed here, that the first saddle-point equation 
Z = Z"^ is obeyed which reduces the number of inte- 
gration variables by a factor 1/2. Configurations which 
are orthogonal to this condition give the same kind of 
factors but have different masses. Doing the Gaussian 
integral js^l we obtain 



In the previous section we have shown that the reduc- 
tion of the theory to a zero mode integral obtains GOE 
spectral correlations. However, we have not yet shown 
under which conditions this reduction is actually legiti- 
mate. This is the question to which we turn next. 



; 7r(m-l)(a+js) i2[ 

mB \ \ 



1 



; 7r(m-l)(s-js) i2 
' mB i 



M I ■ 7r(m-l)(a+jA) l2M 1 „^ 7r(m-l)(s-jA) l2 ' 
L mB -I L 7nB J 



(61) 



where ja = j+ — j- and js = j+ + J- • Differentiating 
w.r.t. the sources we finally obtain the quadratic approx- 
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imation to the correlation function, 

32 



R 



(2) 



1 



E 



(to - l)2(m2B2 



li±=o 

7r2(m-l)2s2) 



(62) 



(to2S2 + 7r2(TO - l)2s2)2 



The contribution of the zero mode (to = 0) is given by 
— -^^Y^ and coincides with the diagonal approximation to 
the GOE correlation function. (Later on we shall see that 
in the case of broken time reversal invariance, one half of 
the matrix elements of Zq become massive implying that 
the contribution of the zero mode reduces to the GUE 
expression 



In the limit B ^ 00, the s-dependence of the contribu- 
tion of massive modes to the correlation function is neg- 
ligible for our purpose, i.e. individual modes contribute 
maximally as '--^ (to— l)2/2m2_B2 ^ [raB)^^ . Only modes 
of mass TO ~ _B~", where a is a non- vanishing positive 
exponent, can survive the limit B 00. The contribu- 
tion of an individual mode is negligible if the exponent 
< a < 1. There are at most 0{B) nearly massless 
modes, and we are led to require that B^"~^ must vanish 
in the limit of large graphs B 00, or that < a < 1/2. 

After these general remarks, let us discuss the masses 
that actually appear in the quadratic action H59(l . We 
first show that modes violating the first saddle-point 
equation Z = Z'^ can safely be neglected. This 
is seen by rewriting the quadratic action as 5(2) = 



TStr 



{Z - Z^)(Z - Z^) + z^z - s+z''siz 

pression shows that fluctuations away from the condi- 
tion Z = Z'^ are suppressed by a large mass of 0(1). 
These fluctuations may safely be ignored, i.e. we may 
assume the condition Z = Z'^ to be rigidly imposed. The 
quadratic action then assumes the reduced form 



This 



5(2) (Z) 



1 



-str 



zz-s+zsiz 



(63) 



where the condition Z — Z"^ reduces the number of inde- 
pendent integration variables by a factor one half. 

We next show that fluctuations Zd,d',d 7^ d' off- 
diagonal in the directional indices may safely be dis- 
carded, too. To this end, let us separate the contribu- 
tion of diagonal and off-diagonal fields, Z'^^^^ and Z°^ , 
respectively, to the quadratic action (|63|l : 



stT{SZS^Z) = str (Z'l'^s Z°«) 



T Q 



(^!?) • (64) 



The matrices Q and 7i contain elements of the type 



Sb^d.b'd'S^^j,,^, or Sb'- 



-d'M^bd^b'd' 



If SLb'd' does not 



vanish there must be a vertex v in the graph, such that 
the directed bond (5', d') ends at v and (5, d) starts at v. 
The partner factors Sb~d.b'd' and Sb'-d'.M then vanish 
(unless the bond 6 is a loop such that (6, d) and (6, —d) 
both start and end at the vertex v. However, for simple 
graphs no loops are present and C/ = 7i = 0.) The ma- 
trix /C contains elements of the form Sb-d.v -d'Sl^y^t 



or Sb'd',bdSbdb'd'- ^'^^ ^ 7^ these vanish (unless the 
bonds b and b' connect the same pair of vertices which, 
however, is forbidden for simple graphs.) For b — b', the 
non-vanishing of the matrix element would again require 
the existence of loops. We thus conclude that /C = 0. 
Decoupled from the scattering operators, the integration 
over modes Zd^d' merely produces a factor of unity (su- 
persymmetry!) so that we will concentrate on the com- 
plementary set of modes 



Zbdd' = Sdd'Zbd 

Zbdd' = Sdd'Zbd, 



(65) 



throughout. The contribution of these conflgurations 
to the generating function is determined by the ele- 
ments of the matrix — {\Sbd,b'd'\'^}- (Here, we used 
that for a time reversal invariant graph, \Sbd,b'd'\'^ — 
Sbd,b'd'S^, i.e. that the matrix T is isotropic in 

time reversal space.) Specifically, the action Sq assumes 
the form 



1 



S(2) = -stT{Z{l-J^)Z) 



(66) 



Within the context of the scmiclassical analysis of ap- 
pendix we have seen that the matrix !F determines 
the classical propagator (the Frobenius-Perron operator) 
on the graph. Comparing with our discussion above, we 
conclude that the eigenvalues of that operator, Ai, deter- 
mine the 'mass spectrum' {rui = 1 — A^} of the theory. 
We have seen that large graphs behave universal if the 
masses scale as rrii ~ B~°', < a < 1/2. Specifically, 
this condition requires the gap 1 — Ai between the zeroth 
Perron-Frobcnius eigenvalue Aq = 1 (corresponding to 
the fully equilibrated zero-mode configuration) and the 
first 'excited' state to scale as 1 — Ai = Ag > i?^",0 < 
a < 1/2. This condition is stricter then Tanners conjec- 
ture < a < 1: For a > 1, corrections to the universal 
result remain sizeable no matter how large the graph is. 
In the intermediate region l/2<a<l — permissible 
by Tanner's criterion — non-universal corrections van- 
ish only if the number r of classical modes with a small 
mass remains constant (or does not grow too fast) such 
that B^A^/r -^00. If, however, the number of low en- 
ergy modes is extensive, r ^ B, the stricter condition 
< a < 1/2 has to be imposed to stabilize universality. 

Above we have shown that in the limit _B — > 00 
only the zero mode effectively contributes to the cor- 
relation function (provided, of course, the master con- 
dition Ag ~ B~" is met.) While the zero mode inte- 
gral must be performed rigorously, all other modes are 
strongly overdamped and may be treated in a quadratic 
approximation. (This is the a posteriori justification for 
the quadratic approximation on which our analysis of the 
mass spectrum was based.) 
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C. GOE-GUE crossover 

The analysis above applied to time reversal invariant 
graphs. In this section we discuss what happens if time 
reversal invariance gets gradually broken, e.g. by ap- 
plication of an external magnetic field. We assume full 
universality, i.e. BA^ ^ 1 such that only the zero-mode 
contributes to i?2(s)- Our aim is to derive a condition 
for the crossover between GOE-statistics (time reversal 
invariance) and GUE-statistics (lack of time reversal in- 
variance.) 

The substructure of the Z-fields in time reversal space 
is given by 



Zdm Zc.h 

yT ^BF jT 



D.b 



7 _BF yT ^ 



7"^ 



(67) 



as a generalization of Eq. (|^ . Here, !Fr> — {\Sbd,b'd'\^} 
while J-c — SM,b'd'S^, (,_^. For a time reversal non- 
invariant graph JFd ^ Tc, and the symmetry of the action 
in time reversal invariance space gets lost. 



Noting that 2B = J2 



bd.b'd' 



\Sbd,b'd'\'^ , we conclude that 
the Cooperon zero mode Zc,b = Yc acquires a mass term 
Bmc str(yc^c)j where the coefhcient 



mc 



2B 



Sbd,b'd'{Sbd,b'd' ~ ^b'd'-^,bd^^' 



bd,b'd' 



(69) 



B 



\tr S^S^af'S^af')\ 



measures the degree of the breaking of the symmetry 
S = . The Cooperon mode may be neglected once 
Bmc — > oo as _B — > oo. 

For the sake of definiteness, let us discuss two concrete 
mechanisms of symmetry breaking: i. breaking the time- 
reversal symmetry of vertex scattering matrices, and ii. 
application a magnetic field. 

Beginning with i., let us consider a large complete 
graph for simplicity. (A graph is complete if any two 
of its vertices are connected by a bond.) The num- 
ber of vertices of these graphs is order B^^"^, and each 



column in Sb'd'M has B^^"^ non-vanishing entries of or- 
der i?^^/'*. Breaking time-reversal symmetry at a single 
vertex thus results in a coefficient mc ^ B^^. Break- 
ing time-reversal invariance at a single vertex is, thus, 
not sufficient to drive the crossover to GUE statistics. 
Rather, a finite fraction ^ B^ {(3 > 0) of time-reversal 
non-invariant vertices is required. Obtained for the sim- 
ple case of complete graphs, it is evident that this con- 
clusion generalizes to generic graphs. 

Turning to ii., the application of a constant magnetic 
field Ab — A causes a global change of all its bond scat- 
tering matrices; We have to replace T{k) — e}-^'^!'^ 



T{k)T{arA) 



A)C/2 ^ This is equivalent to re- 



placing S S{A) = r(crf M)S'T(CTf'M) in the quan- 
tum map. Assuming time-reversal invariance ai A — 
the mass of the cooperon mode becomes 



where Z-^/cm and Z-£,/c.b are 2x2 supermatrices subject mc{A) 

to the constraint Zbb — ^^'^ -^ff ~ ~-^ff' while the 
non-commuting entries of these matrices are independent 
integration variables. The subscripts D(C) allude to the 
fact that in disordered fermion systems, the modes Zd 
(Zc) generate the so-called diffuson (Cooperon) excita- 
tions. Physically, the former (latter) describe the inter- 
ference of two states as they propagate along the same 
path (the same path yet in opposite direction) in config- 
uration space; Cooperon modes are susceptible to time 
reversal invariant breaking perturbations. 

Substituting this representation into the quadratic ac- 
tion, we obtain 

S(') = str (Znil - :Fn)ZD + Zc{l - J'c)Zc) (68) 



1 

2B 
A^H2 



trS't[S'-T(-(T. 
0{A^) 



dir 

3 



A)ST{- 



^dir 



^)]| 



(70) 



where we used aiu^ai — — (T3. We may estimate 



^'■2 



AB 



|tr (^Vf ""^af - 1) I 



(71) 



by setting T{—a^^^ A) k, e ^^2'°'^'" . For a generic scatter- 
ing matrix one may expect |tr {S'^af^'^ Saf^'' — 1) \ ^ B 
such that a small magnetic field of order A ^ B~^^^ 
is strong enough to induce the crossover to GUE statis- 
tics. (Assuming that the geometric 'area' of the graph, 
5 ~ -B, is proportional to the number of bonds, we con- 
clude that the crossover takes place once a finite number 
^ AS xe/h ^ B^ oi fiux quanta pierces the system. This 
crossover criterion is known to apply quite generically in 
disordered or chaotic quantum systems.) 



V. CONCLUSION 

To summarize, we have shown that the two-point spec- 
tral correlation function of individual quantum graphs 
coincides with the prediction of random-matrix theory. 
Corrections to universality vanish in the limit B ^ 00 
provided the gap in the spectrum of the underlying 
'classical' propagator remains constant, or vanishes as 
|m| = |1 - A2I ~ with < a < 1/2. These results 
were obtained by representing the generating functional 
of the two-point correlation function in terms of a non- 
linear cr-model. Closely resembling the theory of spectral 
correlations in disordered fermion systems, this formal- 
ism obtained a fairly accurate picture of correlations in 
the graph spectrum. Specifically, (i) a perturbative ex- 
pansion of the <T-model for large energies establishes the 
contact with semiclassical approaches to the problem, (ii) 
for low energies a non-perturbative integration over the 
fully phase-space equilibrated zero mode configuration of 
the model obtains spectral correlations as predicted by 
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random matrix theory, and (iii) the analysis of the 'mass 
spectrum' of non-uniform modes yields conditions under 
which universality is to be expected: In the limit of large 
graph size, B — > oo, the first non-vanishing eigenvalue 
of the classical Perron-Frobenius operator on the graph 
must be separated from unity by scale Ag larger than 
const. X < a < 1/2. 

This condition turns out to be met by many prominent 
classes of quantum graphs. Examples include complete 
DFT graphs, or complete Neumann graphs It has 
also been shown that almost all unistochastic matrices, 
i.e. matrices of the type J^ij = \Sij\'^, where S runs over 
the unitary group (or, equivalently the circular unitary 
ensemble CUE) , have a finite gap in the limit of large ma- 
trices "33]. For example, star graphs [s^l with the central 
vertex scattering matrix Sij generically display universal 
spectral statistics^ (Counterexamples such as the Neu- 
mann star graph '35] are not generic in this sense.) 

Cases where the universality condition is violated in- 
clude large graphs with low valency (coordination num- 
ber) of vertices. In such systems Anderson localization 
phenomena may in terfere with the buildup of universal 
correlations 



ma y ir 



APPENDIX A: PERIODIC ORBIT THEORY FOR 
GRAPHS AND TANNER'S CONJECTURE 

For completeness we briefly review some elements of 
the periodic orbit approach to spectral statistics on quan- 
tum graphs. Central to the periodic orbit approach is a 
short-time expansion of the spectral form factor, K{t), 
the Fourier transform of the two-point correlation func- 
tion i?2(s) — K{T)e^'^'^^'^ . For moderate bond length 
fluctuations in quantum graphs this quantity is usually 
replaced by the essentially equivalent quantity 



(Al) 



Here n S N is a discrete time corresponding to t = n2B. 
This discrete version of the form factor is connected to 
the correlations in the eigcnphases of the quantum map. 

The form factor HA1|1 is a double sum over periodic 
orbits 7,7' of lengths n where a periodic orbit on the 
graph is a periodic sequence of directed bonds visited. In 
the diagonal approximation — valid for short times n <ti 
2B — only those pairs of orbits are taken into account 
where 7' = 7, or 7' = 7-^ where 7"^ is the time-reversed 
periodic orbit: 



where 



^bd.b'd' 



(A2) 



(A3) 



is the 'classical' propability to be scattered from the di- 
rected bond {b',d') to {b,d). It can be considered as 



the equivalent of the Frobenius-Perron propagator for 
Hamiltonian flows. Due to unitarity of S the matrix F 
is bistochastic J2bd^bd,b'd' = Y^M^b'd'M = 1 and de- 
scribes a Markov process on the directed bonds of the 
graph. The eigenvalues Ai of bistochastic matrices are 
known to lie in the unit circle |Ai| < 1 with at least one 
eigenvalue unity (here corresponding to equidistribution 
on the directed bonds.) 

Universal spectral statistics is expected for 'chaotic' 
graphs in the limit B ^ 00. According to Eq. (|A2() . the 
necessary condition for universality is given by 



tri^ 



t2B 



1, 



(A4) 



where the scaled time t = n/2B <C 1 is kept constant and 
the limit i? ^ 00 is implied. In this case, K'^^'^^ 2t in 

agreement with the short time expansion K'^'^^{t) 
2t of the RMT-form factor 



K 



GOE 



ir) 



(2-log(2|T|- 

-Miogl^ 



1)) 



for |r| < 1 
for Irl > 1. 



Here, r = tA/27r is time measured in units of the RMT- 
level spacing. 

The universality condition above states that any 
propability distribution on the graph will eventually de- 
cay to equidistribution - a Markov process with this 
property is called 'mixing' which implies ergodicity 
(equality of long time-averages to an an average over the 
equidistribution on bonds). This is very week condition 
on a connected graph: a non-ergodic Markov map on a 
graph implies equivalence to a Markov map on a discon- 
nected graph. However, as observed by Tanner ,23i] the 
condition (|A4p is actually stronger than mixing dynam- 
ics (for an example of a mixing graph with non-universal 
spectral statistics — the Neumann star graph, see 35] ). 
This can be seen by rewriting the 1A4|1 in terms of the 
2B eigenvalues Xi of F. Ordering the eigenvalues in mag- 
nitude such that 1 = Al > IA2I > ... and defining the 
spectral gap 



Agap = l-|A2| 



mixing dynamics merely implies A„ 



(A5) 

> 0, i.e. that 
Al = 1 is the only eigenvalue on the unit circle. However, 
there are examples of graphs whose 'classical' dynamics 
is mixing while the form factor does not start as 2r. To 
understand the origin of this exceptional behavior, notice 
that 



tri^" 



2B 



2B 



(A6) 



i=2 



implying the universality criterion | Y^^=2 \ ~^ ^- With 
IE-f2AP;i < E-^2 < (2i3-l)(l-Ag,p)2^^ If 

Agap remains a finite constant as i? — s- 00 this surely 
vanishes and universality is guaranteed. However if 
Agap = ci?~" the correction to unity vanished only if 
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< a < 1. In all known examples of graphs where the 
spectral statistics is non-universal in spite of ergodic clas- 
sical dynamics this condition is, indeed, violated: a> 1. 
For this reason Tanner conjectured that < a < 1 is a 
sufficient universality criterion in the scaling limit r fixed 
while i? — !■ oo. 



The saddle-point conditions identify a zero- 
mode diagonal in directional, and in spin indices 

dmts,d' m' t' s' ^dd' ^mm'^ts^t' s' {^Ob dmts,d'7n' t' s' — 



^dd^ ^mm^^tsM 



) with Y ^Y^ where Y^ = tqY^ t" 



In 

an explicit way of writing, the time reversal structure 
of the Z-matrices (or, equivalently, the F-matrices) is 
given by 



APPENDIX B: TIME-REVERSAL INVARIANT 
GRAPHS WITH SPIN 



In this appendix we discuss the spectral statistics 
of time-reversal invariant graphs with spin (symmetry 
class All, or symplectic symmetry.) This case has been 
considered in connection with the Dirac equation on 
graphs ■ Following a somewhat different approach, 
we will here break spin rotational invariance by choos- 
ing vertex boundary conditions that couple different spin 
components (yet leave time-reversal invariance intact.) 

A spin degree of freedom is straightforwardly intro- 
duced by adding a spin component m — zt^ to the wave 
function on the bonds. This extension turns the quantum 
evolution map Usik) = T{k)ST{k) into a, AB x AB ma- 
trix. We consider non-magnetic graphs {Ab = 0) and as- 
sume independent propagation of the spin components on 
the bonds: T{k)bdm,b'd'm' = Sb,b'Sd,d'Sm,m'e'''^''/'^ . Mix- 
ing of spins occurs at the vertex scattering centers. Com- 
prising all vertex scattering matrices into a single unitary 
matrix S (defined in analogy to spinless case discussed in 
the text) we obtain the condition 



C cT .^dir ^spin rtT spin dir 



(Bl) 



for time-reversal invariance where <J^^^ ^ * ( 1 ) ^^'^ 
Pauli matrix in spin indices m. The corresponding anti- 
unitary time-reversal operator obeys = — 1 charac- 
teristic of symmetry class is All. 

Any system in class All has a doubly degenerate spec- 
trum due to Kramers' degeneracy. By convention, each of 
the doubly degenerate eigenvalues is counted only once, 
i.e. the mean level spacing A = tt/BL as with spinless 
graphs. The oscillatory contribution to the density of 
states is given by Sd{k) — — ^Im ^ln^(fc+). 

In essence, the derivation of field theory representa- 
tion of the generating function parallels the spinless case. 
However, the generalized transposition is now defined as 



Z = 
Z = 



Zc 



Zc 



~-rzi 

z^ 



(B3) 



Projected onto the zero-mode sector, the action is given 
by 

^rsP/,-> 27r p+YY 2ti p_YY 

S^^^(Y,Y) = +i—stT-^ :^-i—stT— — (B4) 

^ ' ' A 1 - YY A 1 - YY 

Except for a factor 2 (which accounts for Kramers' degen- 
eracy) this expression equals the GOE action (|^ : the 
difference between the two cases is hidden in the symme- 
try Y — Y'^. As with the GOE case above, an integration 
over the matrices Y obtains the correlation function 
of the GSE. 

Sufhcient conditions for the zero-mode reducibility 
can be derived as in the spinless case. Modes violat- 
ing the symmetry condition Z = Z'^ , and modes off- 
diagonal in direction space may be discarded. Similarly, 
modes off-diagonal in the spin index m are massive and 
may be neglected, too. Therefore, the validity of the 
saddle-point reduction again relies on the discreteness 
of the eigenmode spectrum of the 'classical' propagator 

^bd 



\m.b' d' m' 



^bdm^b'd'm' 



APPENDIX C: COLOR-FLAVOR 
TRANSFORMATION 

In this appendix we sketch the main conceptual input 
entering the proof of the color-flavor transformation 141(1 . 
(For a detailed exposure of the proof, we refer to .26.]. ) 
Central to our discussion will be a (super)-algebra of 
operators defined as 



[CA,C^] = cac'a - (-) 



_ (_\\A\\A'\ 



CA'CA = Sa,_ 



(CI) 



dir spin 

1 u 



To 



To = - iE^aY + EfaY 



(B2) 



such that T"^ = and = — Tq = f^^- The definition 
of To refiects the different type of time-reversal symmetry. 
Second, the presence of a spin component implies that 
the matrices Zb and Zb — introduced by color flavor 
transformation as before — now have dimension 16 x 16. 
These differences understood, (|47ll applies to graphs with 
spin. 



where A = (a, a), the index a, a' keeps track of all flavor 
components of the theory (boson/fermion, time-reversal, 
directional, etc.), and a, a' = ± distinguish between re- 
tarded and advanced indices. These operators act in an 
auxiliary Fock-space, whose vacuum state is defined by 
the condition 



c+^a|0) = c_^,|0) = 0, 



= 0. 



We may, thus, think of the vacuum as a configuration 
where all (-l-)-states are empty, while all (— )-states are 
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filled; Excitations are formed by creating H — particles and 
(— )-holes. 

The above Fock space contains a sub-space defined 
by the condition that each state contains as many ex- 
cited ( -I- ) -particles as (— )-holes. We call this space as 
the flavor-space. (Alternatively, we may characterize 
the flavor-space as the space of all excitations that can 
be reached from the vacuum state without changing the 
number of particles.) In essence, color-fiavor transfor- 
mation amounts to the construction of two different rep- 
resentations of a projector onto the flavor space. 

Representation no. 1 is constructed as follows: a 
generic state in Fock space can be represented as a linear 
combination |*) = i_ where |^'„+,n_) 

contains particle/hole states. A projection onto 

the flavor state can now be trivially effected by mapping 

Ztt J 

To construct representation no. 2 some more of 
preparatory work is necessary: Consider the Lie- 
supergroup Gl(2i^, 2F) (the group of (4F) -dimensional 
invertible supermatrices.) This group acts in Fock 
space by the representation G1(2F, 2F) B g i-^ Tg = 
exp(cy!i hi{g)A,A'CA'), where Tg acts as linear map in 
Fock space. (For the proof that the assignment g ^ Tg 
meets all criteria required of a group representation, see 

One may show that the representation above is irre- 
ducible in flavor space. This implies that the entire space 
may be generated by the action of the group on a fixed 
reference state, the vacuum state, say. The set of 'co- 
herent states' \g) = Tg\Q) will play a crucial role in the 
construction of the flavor space projector. To bring them 
into a maximally simple form, we first note that transfor- 
mations generated by elements h = bdiag(/i+, € H 
of the sub-group H of matrices block-diagonal in ad- 
vanced/retarded space leave the vacuum invariant (up 
to a constant.) Using the defining commutator rela- 
tions (|Cip it is indeed straightforward to verify that 
Th|Sl) = |ri) X sdet(/i_)^^. Thus, transformations Th do 
not change the vacuum in an 'essential' way and it is suf- 
ficient to consider the action of the coset space G/H on 
1$!). A parameterization of individual cosets gH ^ G/H 
optimally adjusted to our application below reads as 

_/(i-zz)-i/2 z(i-zz)-i/2\ _ 
^ ^ \z{i - zzy^i^ (1 - zz)-i/2 ) - 

_ fi z\ /(l-ZZ)+i/2 ° 

" V V {i-zz)-^/y \z 1 

Letting the transformation Tg act on the vacuum, we 
note that the rightmost factor acts as the identity trans- 
formation while the block-diagonal matrix in the mid- 
dle produces a factor sdet(l — ZZy^^. Finally, the 



left factor is represented by the linear transformation 
exp(c+.QZctQ'C_.ct'), implying that 

\g) = exp(c+,aZ„„,c_,„')l^^> x sdet(l - ZZ)^/^, 
In a similar manner we obtain the conjugate action 

{g-^\ = Sdct(l - ZZ)l/2 X {n\ eM-C^.aZaa'C+,a'), 

where we noted that g~^ — 9\(^z z)^{-z -zy 
claim that 

V"^( dg\g){g\^ 
JG/H 

= j d^{Z, Z)e^>-^°°'=— '|l])(l]|e-^"— ^-'^+-', 

where d^{Z, Z) = d{Z, Z) sdet{l — ZZ), and the flat mea- 
sure d{Z, Z)is the invariant measure on G/K [] is another 
representation of the projector onto the flavor space. To 
prove this statement, we first note that \g) is an element 
of the flavor space (is orthogonal to all non-flavor space 
components contributing to a general Fock space state.) 
Second, V^^ commutes with all transformations Tg: 

T,V"^ I dg'Tg\g'){g'-^\ = 

JG/H 

= 1 dg'\gg'){g'-'\^ j dg' \gg'){{ggr'g\ = 
Jg/h jg/h 

= / dg'\g'){g'-^\Tg^V"Tg. 
Jg/h 

Since the representation g Tg acts irreducibly in flavor 
space, Schur's lemma implies that V^^ must be propor- 
tional to the unit matrix in that space. Finally, the unit 
normalization of V^^ can be shown by computing the 
overlap {n\V^^\n) = J d^{Z,Z) x 1 = 1, where the last 
equality follows from the supersymmetry of the integrand 

We now have everything in store to prove H41|l . Noting 
that = e-xp{c+aV+a+v-aC-a)\^) and (ry+,77_| = 

exp(?7+c(C+Q, -I- C-aV-a) are coherent states, i.e. 

C+aW+, >^-) = +l'+a\v+, V-)-, 

('y+^jy-lc+a = I (+??+«), 

(?7+,77_|c_„ = (77+,r?_| (-r7_„) 

the color flavor-transformation may be proven by map- 
ping the l.h.s. of Eq. (|41|l to a Fock-space matrix element 
and using the equality of the two projector representa- 
tions derived above: 
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J Stt J 2tt 

= {v+,-V-\V'\iy+,iy-) = {v+,-V-\V"W+,iy-) = 

I 
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